Giant Vortices Below the Surface of NbSe 2 Detected Using Low Energy /3-NMR 
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A low energy radioactive beam of polarized 8 Li has been used to observe the vortex lattice near the 
surface of superconducting NbSe2- The inhomogeneous magnetic field distribution associated with 
the vortex lattice was measured using depth-resolved /3-detected NMR. Below T c one observes the 
characteristic lineshape for a triangular vortex lattice which depends on the magnetic penetration 
depth and vortex core radius. The size of the vortex core varies strongly with magnetic field. In 
particular in a low field of 10.8 mT the core radius is much larger than the coherence length. The 
possible origin of these giant vortices is discussed. 



PACS numbers: 74.25.Qt, 74.25.Ha, 76.75.+i, 76.60.-k 

The vortex state of a superconductor exhibits many 
fascinating properties. A fundamental feature, which fol- 
lows from the topology of the superconducting ground 
state and electron pairing, is that each vortex carries an 
elementary quantum of magnetic flux O — hc/2e. In the 
simple Ginzburg-Landau model the order parameter for 
superconductivity is zero at the vortex core and rises to 
the bulk value on the scale of the coherence length, while 
the magnetic field falls away from the core on the scale of 
the London penetration depth. However, even in conven- 
tional superconductors vortices are complex objects with 
properties that were not anticipated when they were first 
predicted by Abrikosov For example, Caroli et al. 
proposed the existence of bound quasi-particles associ- 
ated with an isolated vortex [2j which were eventually 
observed with STM measurement on NbSe2 [![. Ther- 
mal deoccupation of these bound states leads to shrinking 
of the vortex core at low temperatures or Kramer-Pesch 
effect 0. Muon spin rotation (/xSR) has been used ex- 
tensively to measure the magnetic field distribution in 
the vortex state which is a sensitive probe of the vor- 
tex properties. For example /iSR results on NbSe2 [H 
have shown that the low temperature core radius is sig- 
nificantly larger than predicted in the quantum limit [6j . 
Recently there has been considerable work on the role 
of dclocalizcd quasiparticles and the interaction between 
vortices, particularly in multiband superconductors such 
as NbSe2 and MgB2, where there is more than one super- 
conducting gap. It is believed such effects are responsible 
for the unusually large and field dependent specific heat 
[3] and thermal conductivity [8[ in the vortex state. 

In this letter we report observation of giant vortices 
near the surface of superconducting NbSe2 using a novel 
method of low energy /3-detected NMR (/3-NMR) HO. 
The technique is similar in principle to /iSR and in par- 
ticular low energy ^SK [llj , but provides complementary 
information [i"2| . Below T c , we observe a broad asymmet- 
ric /3-NMR lineshape, which is characteristic of a triangu- 



lar lattice of magnetic vortices. As in the case of /iSR the 
internal field distribution, as reflected by the lineshape, 
is sensitive to the magnetic penetration depth and vortex 
core radius. Surprisingly the fitted core radius of 77(10) 
nm is much larger than the coherence length (~ 10 nm) 
in a small magnetic field of 10.8 mT. We propose the ex- 
tended nature of the vortices originates from multiband 
effects and thermal vibrations of the vortices. 

The low energy (30 keV) beam of 8 Li is produced at the 
isotope separator and accelerator (ISAC) at TRIUMF. A 
large nuclear polarization (70%) is generated inflight us- 
ing a collinear optical pumping method. The /3-NMR 
spectrometer sits on a high voltage platform so that the 
implantation energy can be varied between 1 — 30 keV, 
corresponding to an average implantation depth between 
5 — 150 nm. In /3-NMR the nuclear polarization is mon- 
itored through the anisotropic (3 decay. In the case of 
8 Li, which has a mean lifetime of 1.2 s, the emitted be- 
tas have an average energy of about 6 MeV, so that they 
easily pass through stainless steel windows in the ultra 
high vacuum (UHV) chamber. The magnetic resonance 
is detected by monitoring the time-averaged nuclear po- 
larization as a function of a small perpendicular radio 
frequency magnetic field. 8 Li is a spin 1 = 2 nucleus 
with a small electric quadrupole moment Q = +33 mB 
and gyromagnetic ratio 7 = 6.301 MHz/T. In the absence 
of a quadrupolar splitting the lineshape is determined by 
the distribution of local magnetic fields at the Li site. 

The single crystal of 2H-NbSe2 , measuring about 4 mm 
in diameter and 0.1 mm thick, was attached to a sapphire 
plate and mounted on a cold finger cryostat. It had a 
sharp superconducting transition at T c = 7.0 K with a 
width of 0.1 K. The sample was cleaved just prior to in- 
troduction into the UHV (10~ 9 torr); such short exposure 
to air does has no effect on our /3-NMR measurements. 
The beam was focused onto the sample so that there was 
no detectable background signal from the sapphire. 

Fig.QJi shows the /3-NMR spectrum in the normal state 
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of NbSe2 in a magnetic field of 300 mT applied along the 
c-axis and beam direction, which are both perpendicular 
to the surface. Above T c = 7.0 K the lineshape is nearly 
independent of magnetic field, temperature and implan- 
tation depth. In a previous study we showed that the 
Korringa relaxation in the normal state is anomalously 
small compared to simple metals such as Ag 13]. Also, 
there is no resolved quadrupolar splitting. The observed 
line width is attributed to nuclear dipolar broadening 
from the 93 Nb nuclear moments. These results suggest 
that Li occupies a site in the van der Waals gap between 
the NbSe2 layers and is only weakly coupled to the con- 
duction band. 

After field-cooling below T c , the resonance broadens 
asymmetrically (see FigQ})) and exhibits the features 
characteristic of a triangular vortex lattice in a supercon- 
ductor. In particular note the most probable frequency, 
or cusp frequency, shifts by an amount A c below the 
normal state frequency. The cusp frequency (or magnetic 
field) arises from Li located midway between two adjacent 
vortices. In addition, there is a high frequency tail which 
corresponds to the magnetic field distribution near the 
vortex core. The high frequency cutoff, which is shifted 
by an amount A„ above the normal state frequency, cor- 
responds to the field at the vortex core. At high implan- 
tation energies we expect the magnetic field distribution 
(frequency distribution) to be close to that in a bulk su- 
perconductor. Such bulk field distributions have been 
studied extensively with muon spin rotation (/iSR ) and 
used to extract the magnetic penetration depth and vor- 
tex core radius. In applied magnetic fields (H << H c2 ) 
the local magnetic field at position r = (x, y) relative to 
a vortex can be decomposed into its Fourier components: 
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(1) 



where the sum is over all reciprocal lattice vectors of the 
vortex lattice, A 6 is the in-plane penetration depth, p is 
the vortex core radius, Bq is the average magnetic field, 
and C(k,p) is a phenomenological cutoff function which 
characterizes the shape and size of the vortex core[LJ|. p 
is a function of the coherence length, but may also de- 
pend on the electronic structure of the vortex and vortex- 
vortex interactions. The best fit[l5[ of the current data 
is with a simple Gaussian cutoff C(k,p) — exp [— \k 2 p 2 ] 
which gives A a {,=230(30) nm and p=13(l) nm. These 
parameters depend slightly on the model for C(k,p). 
For example using a Bessel function cutoff 16] gives 
X a b = 230(30) nm and p = 11(1) nm. Previous /iSR work 
on bulk NbSe2 [l9| obtained X a b ~ 170 and p w 12 at the 
same field and temperature, using a Bessel function cut- 
off. The agreement is reasonable considering there are 
substantial differences in the experimental methods, the 
form of the raw data and resulting analysis. For exam- 
ple, /J-NMR spectra are acquired in the frequency domain 
rather than in the time domain for /iSR . Thus impor- 
tant features such as and A c are evident in the raw 
/3-NMR data without data processing or analysis. Also, 



in /3-NMR there is no detectable background that inter- 
feres with the signal. On the other hand the form of the 
field distribution is more complicated as explained below. 
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FIG. 1: (Color online) (a)The /3-NMR spectrum in the normal 
state of NbSe2 at 10 K in a magnetic field of 300 mT applied 
along the c-axis. The beam energy of 20 keV corresponds 
to a mean implantation depth (z) — 85 nm. The observed 
Gaussian line broadening is attributed to 93 Nb nuclear mo- 
ments, (b) The same conditions as (a) except field cooled to 
3.5 K or 0.5T C . The asymmetric lineshape is characteristic of 
a triangular lattice of magnetic vortices. 

Fig. [2] shows /3-NMR spectra measured in a lower 
magnetic field of 10.8 mT. The normal state resonance 
(Fig. [2^) is very similar to that observed at 300 mT. The 
spectra in the superconducting state (Figs. [5Jd and [2fc) 
show a strong dependence on temperature (not shown) 
and implantation depth (z). Several important features 
are evident from Figs. [2}} and [2J: without any fitting. 
Firstly, A c decreases as the mean depth changes from 
85 nm (Fig. [2t) to 8 nm (Fig. [2b). Such narrowing of 
the field distribution near the surface is understandable 
since the vortex lattice must approach that of a bulk 
superconductor deep inside the material; whereas, any 
line broadening from the vortex lattice must vanish well 
outside the material. The actual crossover occurs when 
the mean depth becomes comparable to ao/2n where 
ao[nm] = 1 546/ \J B ext [mT] is the mean spacing between 
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vortices [17j. Secondly, A c in Fig. [2J: has increased by 
about 50% compared to the corresponding lineshape at 
300 mT shown in Fig.[]i>. This is also consistent with the 
vortex lattice model since the vortices are further apart 
in low field. 
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FIG. 2: (Color online) /3-NMR spectra in NbSe2 in a low mag- 
netic field of 10.84 mT. (a) In the normal state at 10 K with 
a beam energy which corresponds to a mean implantation of 
(z) = 8 nm. (b) The same conditions as (a) except field cooled 
to 0.5T C . (c) The same temperature and magnetic field as (b) 
but the mean implantation depth is about 10 times larger. 

The most unexpected feature in the data is the re- 
markable similarity between the lineshapes in Figs. [TJd 
and [2}}, given the magnetic fields differ by a factor of 27. 
In particular, in Fig. [5b note the cutoff at a frequency 
A„ above the normal state frequency. Recall, this cor- 
responds to the frequency (magnetic field) at the vortex 
core. For example using p = 13 nm and X a b = 230 nm 
one finds A„ = 76 kHz or about three times larger than 
what is observed. Typically A„ can only be observed 
in much higher magnetic fields such as in Fig. [TJd where 
the vortices are close together. The small value of A„ 
in Figs. and[5fc indicates a very extended vortex core 
where p 3> £ a fc- The possible origin of the over-sized 
vortices is discussed below. 

All the data in the superconducting state were fit to a 
vortex lineshape model, which includes effects from the 



surface. Following Ref.[l7|, we assume Laplace's equa- 
tion is valid outside the superconductor, and that a mod- 
ified London model for a triangular lattice of vortices ap- 
plies inside the superconductor. This yields a magnetic 
field at a depth z equal to: 



B(r,z) = B a J2 
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where A 2 = k 2 + 1/A 2 fc . In analogy with bulk super- 
conductors we include a simple Gaussian cutoff function 
C(k,p) — exp [— \k 2 p 2 \. For each energy we calculate a 



depth-averaged field distribution: 

(n(S)> - / f(z)n(B,z)dz. 



(3) 



where f{z) is the stopping distribution obtained using 
the TRIM.SP code to simulate the implantation of 8 Li 
in NbSe2 [H| and n(B, z) is the field distribution at a well 
defined depth z. The final step is to convolute the result- 
ing frequency spectrum with a Gaussian broadening func- 
tion which takes into account the nuclear dipolar fields 
and any residual disorder in the vortex lattice [iH]. Typi- 
cal fits to the model field distribution are shown in Fig. [2] 
The average frequency equals the normal state frequency 
implying there is no measurable flux expulsion, which is 
reasonable for our sample geometry. The measured val- 
ues for A c and A„ are shown in Fig. [3l For comparison 
the solid lines are the best fit to the vortex lineshape 
model assuming a depth independent X a b and p. At 10.8 
mT we obtain X a b — 167(15) nm and p = 77(10) nm 
compared to A a & = 230(30) nm and p = 13(1) nm at 300 
mT. At both fields A c and A„ behave roughly according 
to this simple model, but there are clear deviations. In 
particular A c increases with decreasing (z) close to the 
surface, whereas the simple model predicts a monotonic 
decrease. Also, both A c and A„ show a stronger variation 
with depth than the simple model predicts. These are in- 
dications that the vortex interactions depend slightly on 
depth. 

The fitted value for X ab at 10.8 mT is about 30% 
smaller than at 300 mT and is very close to the /xSR 
results for bulk NbSe2. However, the most significant ef- 
fect of the magnetic field is the extremely large value of p 
at 10.8mT. At 300 mT p is 13(1) nm, which is close to the 
expected coherence length and measurements of the core 
radius made with /j,SR in higher fields [l9| . At 10.8mT the 
effective p increases to 77(10) nm. It should be stressed 
that changing the functional form of the lineshape pro- 
duces only small changes in p and has no influence on the 
ratio between p in high and low field. The magnitude of 
the effect is much larger than the calculated field depen- 
dence of the core radius originating from the bound states 
which cause the Kramer-Pesch effect [20| . It is interesting 
to compare these results with /iSR in bulk NbSe2. Early 
studies in this temperature range showed a substantial 
increase in the vortex core size as the field was lowered 
but the signal was not followed below 100 mT [lj|. The 
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FIG. 3: (Color online) Depth dependence of the cusp fre- 
quency A c and vortex core frequency A„ relative to the nor- 
mal state frequency. The open circles are measured in a mag- 
netic field of 10.8 mT while the filled circles are taken at 300 
mT. In all cases the temperature is 3.5K which corresponds to 
0.5 T c . The solid curves are from the vortex lineshape model 
with a single depth independent A a ;, and p. 



field dependence we observe indicates a stronger varia- 
tion with field below 100 mT than predicted from Ref. 
fl9| where the extrapolated low field value for the core 
size was 25 nm at 0.5T C . Also, very recent work at much 
lower temperatures (20 mK) gives an extrapolated low 
field value of 10 nm [2l| . A 60% increase in p below 400 
mT is attributed to a multiband effect, whereby a sec- 
ond smaller superconducting gap determines the vortex 
radius in low fields. In a simple model, with minimal 



band coupling, the low field £ a b w vf/ttAq where vf is 
the Fermi velocity and Ao is the smaller gap. We observe 
a much larger field dependence to the core radius. The 
current measurement was made in a lower magnetic field 
and of course much closer to the surface, but the key dif- 
ference is likely the temperature. Both ARPES [22| and 
STM[23| find evidence for a broad distribution of gaps 
with parts of the Fermi surface showing no measurable 
gap at higher temperatures near T c . The ratio between 
the maximum and minimum gaps is about five at 0.5T C 
[23l ] . Thermal vibrations of the vortex lattice may also 
contribute to the observed field dependence in the core 
radius, since in low magnetic fields the vortices are fur- 
ther apart and less confined by their neighbors. The up- 
turn in A c for small (z) (see Fig. [3]) indicates some depth 
dependence to vortex interactions as one approaches the 
surface. This is consistent with vibrations since the vor- 
tex lattice is expected to stiffen very close to the surface, 
either due to pinning or increased vortex repulsion [24J. 
A comprehensive theory of the vortices and their inter- 
actions, which includes both electronic and vibrational 
excitations is needed to resolve the precise origin of this 
effect. 

In conclusion, low energy /3-NMR has been used to 
measure the field distribution of a vortex lattice below 
the surface of a superconductor. Below T c in NbSe2 the 
resonance line is inhomogeneously broadened and has a 
lineshape which is characteristic of a triangular vortex 
lattice. In a magnetic field of 300 mT the magnetic pen- 
etration depth and core radius are similar to //SR mea- 
surements in bulk NbSe2. However, in a lower field of 
10.8 mT the vortices are much larger. The multiband 
nature of the superconductivity and thermal vibrations 
of the vortices are a possible explanation. 
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